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Nuclei colliding at very high energy create a strong, quasi-classical gluon field during the initial 
phase of their interaction. We present an analytic calculation of the initial space-time evolution 
of this field in the limit of very high energies using a formal recursive solution of the Yang-Mills 
equations. We provide analytic expressions for the initial chromo-electric and chromo-magnetic 
fields and for their energy-momentum tensor. In particular, we discuss event-averaged results for 
energy density and energy flow as well as for longitudinal and transverse pressure of this system. For 
example, we find that the ratio of longitudinal to transverse pressure very early in the system behaves 
as Pl/pt = —[1 — where r is the longitudinal proper time, Q is 

related to the saturation scales Qs of the two nuclei, and a = hi{Q^/m^) with m a scale to be defined 
later. Our results are generally applicable if r < 1/Q. As already discussed in a previous paper, 
the transverse energy flow S* of the gluon field exhibits hydrodynamic-like contributions that follow 
transverse gradients of the energy density VY. In addition, a rapidity-odd energy flow also emerges 
from the non-abelian analog of Gauss’ Law and generates non-vanishing angular momentum of the 
field. We will discuss the space-time picture that emerges from our analysis and its implications for 
observables in heavy ion collisions. 


I. INTRODUCTION 

Collisions of nuclei at high energy at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Col¬ 
lider (LHC) have established the existence of a deconfined 
phase of partons at high energy densities e > 1 GeV/fm^ 
m- The future goal of these programs is to make pre¬ 
cision measurements of properties of quark gluon plasma 
(QGP) and to study further details of the phase dia¬ 
gram of quantum chromodynamics (QCD). This ambi¬ 
tious task requires a detailed understanding of the bulk 
dynamics in nuclear collisions. The most promising can¬ 
didate theory for understanding the initial phase of these 
collisions is color glass condensate (CGC) [IHZ] in which 
the initial interaction of nuclei, modeled as a collection 
of SU{3) color charges before the collision, leads to a 
quasi-classical gluon field after the collision. This field 
eventually decays into a thermalized QGP. 

Once the system is close to local kinetic equilibrium, 
dissipative relativistic fluid dynamics has become the tool 
of choice to compute the expansion and cooling of the 
QGP fireball. Comparisons of hydrodynamic simulations 
with experimental data have been increasingly successful 
in pinning down the shear viscosity and the equation of 
state of high temperature nuclear matter [SHE]. The 
equilibration time Tth, when hydrodynamic concepts can 


be applied, as well as the initial values for energy den¬ 
sity, energy flow, and all other currents at r = Tth, are 
often treated as parameters in the fluid dynamic simula¬ 
tion. Model calculations of the initial state, such as the 
Glauber model [T^, often constrain only a small subset 
of initial parameters. In particular, initial transverse flow 
is still often poorly constrained in many calculations or 
even neglected despite very good arguments to the con¬ 
trary [Hill]. If color glass condensate is found to be 
the applicable description of the initial interaction of nu¬ 
clei at collider energies we will, in principle, be able to 
calculate the initial conditions at the time Tth. Recent 
progress seems to indicate that this is the correct path 

[Mini- 

Here we have a modest goal. We would like to present 
analytic results that bridge the gap between known re¬ 
sults for the classical gluon field in single nuclei before the 
collision [TS] and the glasma helds at a time tq ~ 1/Q 
after the collision. The tq represents the limit of conver¬ 
gence of the small-time expansion we employ. However, 
in terms of physics it also represents the point at which 
the longitudinal pressure pp, initially large and negative, 
approaches zero or even becomes positive, a necessary 
(but not sufficient) condition for pressure isotropization. 
Our results then provide solid and urgently needed input 
to constrain the energy-momentum tensor at a later time 
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Tth > tq which can feed into fluid dynamic simulations. 
It might be used in an ad-hoc thermalization approxima¬ 
tion, as in HZIIIS], or it might serve as the starting point 
of further studies of thermalization itself HOI [SB- The 
phenomenology we find is surprisingly rich. For example, 
the system has non-zero angular momentum and exhibits 
directed flow. It resembles aspects of phenomenological 
models based on QCD strings or string ropes suggested 
previously |221123] ; however, our derivation here is based 
strictly on classical QCD. 

We should note that while tq is rather early in the col¬ 
lision, it is within this initial time period that important 
global properties are set. These include how much en¬ 
ergy, momentum, and angular momentum are transfered 
from the initial system of colliding nuclei and deposited 
in the relevant part of the fireball around midrapidity. 
While we will focus on analytic results for event-averaged 
quantities, it is in principle straight-forward to construct 
a semi-analytic event generator based on our results. 

Color glass condensate has been developed from the 
idea that nuclear wave functions in the asymptotic limit 
of very high energies should exhibit novel properties of 
QCD [THl This state is characterized by a 

slowing growth of the gluon distribution with increasing 
energy (or decreasing Bjorken-x). The gluon area density 
in a hadron or nucleus saturates and thus defines a sat¬ 
uration scale Qs- We will denote the proper saturation 
scale in a nucleus by Qs and will assume it is related to 
the scale Q used earlier by a numerical factor. We will 
discuss the ultraviolet scale Q in more detail later. At 
high energies Qs becomes large, Qs S> Aqcd, and the 
strong coupling as becomes small. The Qs is assumed 
to be on the order of a few GeV in heavy ion collider 
experiments. In addition, gluon occupation numbers are 
large and a quasi-classical description of the gluon field 
becomes applicable. If two nuclei collide at high en¬ 
ergy, the interaction of the two color glass condensate 
states create what is sometimes referred to as glasma 
[n m nil [Ml US- Here we are interested in the early 
time evolution of glasma. We will use the classical ap¬ 
proximation, known as the McLerran-Venugopalan (MV) 
model UlillMlIli]- We will, however, need to general¬ 
ize the original form of the MV model in this work to 
allow for a rigorous description of transverse dynamics. 
Quantum corrections have been studied and seem to indi¬ 
cate that the classical description is adequate to describe 
the evolution of the system up to times of order 1/Qs 
[20] . Initial small fluctuations can grow exponentially at 
times beyond 1/Qs and lead to instabilities. They are 
probably an important step on the path to thermaliza¬ 
tion. Recently, important progress has been made on this 
phase in the evolution of gluon fields [201 El] ■ 

The time tq has multiple important implications in 
our work. It signals the breakdown of the classical ap¬ 
proximation as well as the limit (on purely mathemat¬ 
ical grounds) of our specific solution to the Yang-Mills 
equations. However, it also heralds decoherence of the 
classical fields [30] at which the net transfer of energy 


and angular momentum from the receding nuclei onto the 
fireball presumably stops, and it is responsible for most 
of the reduction of the pressure asymmetry (neglecting 
transverse gradients) 


Pl -Pt 


{pL + 2 pt)/5 


= -6 






+ 0{Qt)\ (1) 


where a = ln{Q^/ tjt). This will be discussed in Sec. VI 


Our paper is organized as follows. In Sec. |TT| we re¬ 
view the MV model for single nuclei and colliding nu¬ 
clei on the light cone. We discuss a recursive solution 
of the equations of motion of the gluon field. We also 
compare the emerging space-time picture to existing phe¬ 
nomenological approaches. In Sec. HI we calculate the 
energy-momentum tensor of the early gluon held as a 
function of the initial color electric and magnetic helds up 
to fourth order in proper time r. In Sec. IV we general¬ 
ize the assumptions used to calculate expectation values 
of observables in the MV model and redo the classical 
calculation of the gluon distribution function of a nu¬ 
cleus. We then proceed to calculate the expectation val¬ 
ues, or event averages, of gluon held correlation functions 
of higher twist which will be needed later on. In Sec. V 
we compute the expectation value of the glasma energy- 
momentum tensor up to fourth order in r, although at 
third and fourth order in r only leading contributions in 
Q are computed. In Sec. VI we explore the phenomeno¬ 
logical consequences including pressure anisotropies and 
how. Section VH summarizes our results. 


II. THE GLUON FIELD OF TWO COLOR 
CHARGES ON THE LIGHT CONE 

In this section we discuss analytic solutions for the 
Yang-Mills equations of two nuclei colliding on the light 
cone with non-abelian SU{Nc) charges kept hxed. The 
setup is reminiscent of an expanding color capacitor: in- 
hnitely Lorentz-contracted sheets of SU{Nc) color charge 
move towards each other (along the z axis), pass through 
each other, and recede. Color capacitor-like systems have 
been discussed in the literature in other contexts and we 
will come back to a comparison later on. The CGC setup 
is briehy reviewed in the following. 

In the CGC limit nuclei move on the light cone. Their 
partons can be divided into source partons with large 
momentum fraction x and classical gluon fields that ef¬ 
fectively describe small-x gluons in the nuclear wave 
functions, as first discussed by McLerran and Venu- 
gopalan [HE]- The source is given by a SU{Nc) cur¬ 
rent Note that we use underlined upper 

or lower indices for SU{Nc). The t- are the Gell-Mann 
matrices. We have specified our definitions in appendix 
[A| The gluon field strength and its gauge held 
couple to the current through the Yang-Mills equations 


( 2 ) 









A. General Shape of the Field 
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and the continuity equation 

[D^,J^]=0. (3) 

The internal dynamics of the source are frozen on time 
scales that describe interactions with probes or other nu¬ 
clei (the glass in CGC) and are therefore kept fixed on 
the light cone. In addition, during a collision large angle 
scatterings of source partons are rare (those would be re¬ 
ferred to as hard processes). The slowing down of source 
partons through the interaction, in other words the back 
reaction of the field on the sources, can be significant but 
at sufficiently large collision energy the source partons are 
close to the light cone even after the collision. This has 
been confirmed experimentally even at top RHIC ener¬ 
gies where nuclei, represented by the net baryon number 
carried by the valence quarks, lose about three quarters 
of their kinetic energy during the collision. Therefore 
they stay ultrarelativistic throughout [3T]. This justifies 
the assumption that a current along the -|- light cone is 
invariant, or independent of the x'^ coordinate. Because 
of the practically infinite Lorentz boost the source is also 
infinitely thin in the x~ direction, and can therefore be 
solely described by an S'C/(3)-valued area density p{x±), 
where x^_ is the vector of transverse coordinates. For our 
definitions of light cone coordinates we refer the reader 
to appendix [5] 

We represent two colliding nuclei on the light cone 
through two currents Jf 2 along the -I- and — light cone, 
respectively, given by two SU{3) charge densities pi{x±) 
and p 2 ix±). The components of the currents in light cone 
coordinates are 

Jiix) = S{x~)pi{x±), Jf(x)=0, (4) 

J 2 ix) = S{x~^)p 2 ix±), J^{x)=0, (5) 

Jl2ix) = 0. ( 6 ) 

with i = 1,2. The total current Jf -|- satisfies the 
equation of continuity ^ if we choose an axial gauge 
with 

x'^A~ +x~A+ =0. (7) 

We will keep this choice of gauge throughout this section. 

We note that nuclei fixed on the light cone will lead to a 
boost-invariant system after the collision. In particular, 
the energy-momentum tensor of the gluon field will be 
boost-invariant. This will be an important caveat when 
we discuss the global space-time structure of the fireball. 
In our calculation global energy, momentum, and angular 
momentum are not conserved as the nuclei are reservoirs 
for those conserved quantities. In reality, those quantities 
are finite and conserved. However, we still expect to 
gain realistic insights of the rapidity densities of those 
quantities as long as we stay far enough away from the 
final rapidities of the nuclei. Corrections to the boost 
invariant approximation can in principle be taken into 
account [5^ . 


Kovner, McLerran and Weigert were the first to discuss 
the general space-time structure of the gluon field in the 
CGC formalism in the collision of two nuclei [23]. One 
can write down the following ansatz for the a;"'"-a:“-plane: 

A'^{x) =0(a;^)0(x“)x"''A(T, x_l) , 

A~ (x) = — 0(x~^)0(x~)x~A(t, x±) , 

A^(x) =0(x~)0(—x~^)Al(x±) (8) 

-I- 0(a:''’)0(—a:“)A 2 (af_L) 

-b 0(a;+)0(a;“)A(L(T,f_L), 


where again i = 1,2. The A\{x±_) and A|(af_L) are 
the gluon fields of the single nuclei before the colli¬ 
sion, which are purely transverse in this gauge. Here 
r = - is the longitudinal proper time. The 

A(t, x_l) and A\{t,x±) are smooth functions in the for¬ 
ward light cone and describe the field after the colli¬ 
sion. They are the glasma fields we will be interested 
in. There is no explicit dependence on the space-time 
rapidity 77 = \ \n{x'^/x~) in A and A\, reflecting the 
boost-invariance of the system. Figure shows the dif¬ 
ferent regions of the light cone including the region of 
applicability of this work. 

In each sector of the light cone the Yang-Mills equa¬ 
tions have to be satisfied separately. In the forward light 
cone they can be written in the convenient form |24j 





\ d 1 

1 

r ■ a 1 

igr 


r 



( 9 ) 

( 10 ) 


1 ^ d .i 

T OT OT 


igr^ [A, [D\A]\ - [D\FA]=Q. 

( 11 ) 


The field strength tensor in the forward light cone can 
be expressed in terms of the gauge potentials A and A\ 
in this gauge as 


A1 T [D\ A]^ , (12) 

= d^A]_ - d^A]_ - ig[A^^,A\}. 

Boundary conditions connect different light cone sectors. 
The ones for the forward light cone read [24] 

A\{t = {),Xi_) = A\{xx) + A\{xi_), ( 13 ) 

A{t = 0 ,£_ l ) = [A\{xi_),A\{xi_)\ . ( 14 ) 

We interpret them as initial conditions for the fields at 
T = 0 for the fields in the forward light cone r > 0. 
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FIG. 1. (color online) The ^-i-plane with the two currents Ji and J 2 given on the x^- and a:“-axes, and the four regions given 
by the solution Q. The theoretical limit of the classical approximation in the forward light cone at r = tq and the approximate 
thermalization time rth are shown schematically. 


Equations (|9) through ( |11[ ) together with the condi¬ 
tions (13) and pose the boundary value problem to 
be solved. An analytic solution in closed form is not 
known for the most general case. The weak field or 
abelian limit was first treated in |24j and will be repro¬ 
duced below. Several groups have discussed numerical 
solutions [roi[55H55] . usually focusing on the plane ry = 0. 


A different approach to solve the problem was hrst ad¬ 
vocated by some of us in [191 ES] ■ The basic idea is as 
follows. Since the classical approach to CGC loses its 
applicability very soon after the collision, it will be suffi¬ 
cient to focus on the near-field, or small proper times t. 
In that case one can utilize a systematic expansion of the 
Yang-Mills equation in a power series in r [50115S] . We 
can expect to find the leading terms in such an expan¬ 
sion analytically. The natural scale for the convergence of 
such series should be given by the only time scale in the 
problem, namely, 1/Qs- We will see that this is indeed 
the case. 


B. r-Expansion and Recursive Solution 

Let us dehne the power series 


x±) = (£_ l ), (15) 

n—0 
00 

= (16) 
n=0 


for the fields parameterizing the gauge potential in the 
forward light cone. We devise equivalent power series for 
the field strength tensor, covariant derivatives and the 
energy-momentum tensor. We do not include any diver¬ 
gent (1 /t") or logarithmic (Inr) terms in r. While the 
field equations themselves can have divergent solutions 
they have to be discarded because of the boundary con¬ 
ditions ([T^ and (jldl). 


We can discuss this point in more detail for the abelian- 
ized version of the equations. In the case of weak helds 
the non-linear terms in the Yang-Mills equations are usu¬ 
ally neglected, leading to a greatly simplified abelian ver¬ 
sion of the boundary value problem. The analytic so¬ 
lution in closed form can be readily found [53]. After 
applying a Fourier transformation of t he transverse co¬ 
ordinate, 9* —)■ Eqs. ^ and (11) take the form of 

Bessel equations 

1 . 1 d 1 1 ,1 n 

—-|- —^ ——zA —zA — —irzA — 0, (1^) 

z az^ dz z z'^ 

z':^ai + z^ai + z2ai = 0, (18) 

dz^ dz 

where z = k±T. A physical polarization V*A(|_ = 0 has 
been chosen for the transverse field. There are two inde¬ 
pendent sets of solutions, Bessel functions of the first kind 
A ^ Ji{z)lz, Af ^ Jo{z) which are regular at r = 0, and 
Neumann functions A ^ 7Vi(z)/z, A]_ ~ Ao(z) which 
lead to singular solutions A ^ z~^, Af ^ Inr for r —>■ 0. 
The solut ion with Neumann functions is not compatible 
with Eq. (10) which imposes djdrAf = 0. The singular 


solution therefore has to be excluded. 
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Let us now return to the solution of the general non- 
abelian problem. The power series turns the set of 3 
differential equations 0. ( [To| ), and ( [TT| ) in x± and r 
into an infinite system of differential equations in x_\_- 
Amusingly, we can solve this system recursively. The 
boundary conditions (13) and (14) provide the starting 
point of the recursion 


Ai — 4* 

^-L(0) “ ^1 + /I 2 , 




(19) 

( 20 ) 


It can be shown that all coefficients of odd powers vanish, 
■^{ 2 k+i) = 0 and A(L( 2 fe+i) = 0- Finally, one finds the 
recursion relations for even n, n > 1, to be 


^(n) / I r)\ 

' ' n{n + 2) 


E 


D) 


(*)’ 




^(l) > ^(m) 


„2 ( E 

\k-\-l—n—2 


T)3 pJ'i 


( 21 ) 


+ ig 


E 

k-\-l-\-7n—n—A 




One can readily see that these expressions solve ([^ and 
(11). It is less straight-forward to show that the recursion 


relation solves Eq. (10). One can go order by order in r. 


and we have explicitly checked that our recursive solution 
solves Eq. (10) up to 4th order in r. 


One can use the abelianized case for a cross check. 
After dropping non-linear terms, and after applying a 
Fourier transformation to the transverse coordinates, the 
recursive solutions can be easily cast in the form 


(">1) (22) 


1 


aLOi _ \n /2 AhOi 

2^-L(n) 7^112 1 ^-L(0) ’ 


(23) 


where the double factorial is nil = n{n — 2)(n — 4) • • • 
and the index LO signals the abelian case. These terms 
are just the coefficients of the Bessel functions already 
discussed above. 


A^°(r,fci) (^kxT) , 

rC_l_ T 

(24) 

A^°*(T,fcj_) =A^‘^J)(kj_)Jo ik±T) . 

(25) 


Thus the small-r expansion immediately recovers the full 
abelian solution. 


The recursive solution (21) and its consequences are 
the basis for the remainder of this manuscript. A brief 
discussion on the convergence of the series expansion is 
in order. From the abelian case above we infer that the 
radius of convergence is 00 , independent of the charge 
distributions /9i,2; as long as we are in the weak field 
limit. In the opposite limit of extremely strong fields one 
can make the following estimate. Keeping only the max¬ 
imally non-abelian terms, we expect from the recursion 
relations that 


lA")! ~ l'724r+Vl and ~ \gA\^\A \, (26) 


where 


|A| = /4iE> (27) 

is written in terms of the fields in the initial nuclei before 
collision. Here we assume head-on collisions of equal nu¬ 
clei (pi = P 2 ) for simplicity of argument. We anticipate 
from our results later on that \gA\'^ ^ g'^p/A'K, where p 
is the color charge density of the two incoming nuclei. 
From the geometric interpretation of the saturation scale 
Qs we further have ^ [35]. Hence we find that 

parametrically 

(28) 

\A^xin)\ - Q^\A\ ■ (29) 

This suggests that the convergence radius of the series 
in this extreme case is indeed parametrically set by the 
saturation scale, Tconv ~ l/Qs- We can find further phe¬ 
nomenological validation in Sec. |VI A| when we compare 
to numerical solutions of the Yang-Mills equations. 


C. The Near Field 


A resummation similar to the abelian case seems elu¬ 
sive for the general solution. However, we can analyze 
the few lowest order terms explicitly. This amounts to a 
description of the near field close to the light cone. The 
series expansions for the gauge potential are 


AiT,xx) = A(o) + ^[D=, A(o)]] + [D\ [D^, A(o)]]]] + ^ 


^[D^A^o),D^Bo]+0{t^), (30) 


Al(r,x^) = + ^e^WW'^D'^Bo - '-^[Bo,D^Bo] + ^[A(o), [D\A^o)]] + 0 (t^} , 


16 


64 


(31) 
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FIG. 2. (color online) Left: After the collision Lorentz-contracted nuclei with color charges and transverse fields develop 
longitudinal fields Eq and Bo between them. Right: Transverse fields between the nuclei are induced by the decrease of 
longitudinal fields after a short time. Only fields from Faraday’s and Ampere’s Law are indicated. 


where we have used the short hand notation = 

d'‘ — In the remainder of this work will denote 

the covariant derivative with respect to the initial gauge 
field and we will mention explicitly if we refer to covariant 
derivatives at other times. The Bq is the longitudinal 
chromo-magnetic field which is discussed below. 

Let us carry out an order by order analysis for the held 
strength tensor 

E =E(o) + tE(i) + t^E(2) + ... , (32) 

E + 'rB(i) + 'r^B(2) + .. ■ , (33) 


of chromo-electric and chromo-magnetic helds. From 
here on electric and magnetic always refer to chromo¬ 
electric and chromo-magnetic. The components of the 
held strength tensor can be readily computed from the 
gauge potential using Eqs. (12). We observe that only 
the longitudinal components of the electric and magnetic 
helds have non-vanishing values at r = 0 (TH] 


Eq = 

Eq = B^o) 


- ^(0) 


= F, 


21 

(0) 


= ig [-41,^2] , 


= igi 


ij 




(34) 

(35) 


They can be seen as the seed helds for the glasma de¬ 
veloping in the forward light cone. The transverse helds 
vanish at r = 0: = 0. 

The dominance of longitudinal helds, both electric and 
magnetic, at early times has been discussed in [28l |29] 
and has since been often relabeled as the occurrence of 
color hux tubes. They are similar but not directly compa¬ 
rable to QCD strings. QCD strings are a reaction of the 


QCD vacuum to color charges. Here we consider helds 
close to the center of a collision of large nuclei which are 
far removed from the QCD vacuum. Non-trivial QCD 
vacuum effects are not included in the classical Yang- 
Mills picture considered here. QCD strings have been 
successfully used to describe collisions of nucleons at large 
energies m- It would be desirable to hnd a natural tran¬ 
sition between glasma hux tubes in the center of collisions 
and QCD strings describing the dynamics at the bound¬ 
ary of the collision zone, but that is beyond the scope of 
this work. The initial longitudinal magnetic and electric 
helds can be of similar strength in the glasma. Figure 
shows a sketch with nuclei consisting of Lorentz con¬ 
tracted sources and transverse gluon helds, and longitu¬ 
dinal helds stretching between them after the collision. 

It is useful to briehy mention a reinterpretation of the 
initial longitudinal helds pointed out in |28j . This can 
help to make connections with some existing phenomeno¬ 
logical models in which an exchange of color charge be¬ 
tween the nuclei is envisioned at the time of their over¬ 
lap [22] [23] . Those effective color-anti-color charges on 
opposite nuclei then lead, in a quasi-abelian picture, to 
longitudinal electric helds between the nuclei after they 
have separated and recede from each other. This ap¬ 
pears similar to the early glasma picture. Note, however, 
that here the charges pi and p 2 are strictly kept constant 
throughout the collision and the longitudinal held arises 
from non-abelian interactions between the helds of the 
two nuclei. But to aid our intuition, we can rewrite the 
covariant derivatives in Gauss’ Laws for chromo-electric 
and chromo-magnetic helds as ordinary derivatives and 
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FIG. 3. (color online) Transverse fields in the 2 ;-j/-plane for simple abelian analogue to Eqs. (371, (38 1 for randomly simulated 
initial fields Al and A 2 . Upper panels: Electric field E‘ (arrows) on background of the longitudinal magnetic field Bo (shading). 
Lower panels: Magnetic field B‘ on background of Bq. Left panels: Rapidity-even terms from Faraday’s and Ampere’s Law, 
respectively (corresponding to ?7 = 0). Middle panels: Rapidity-odd terms from Gauss’ Law. Right panels: Full transverse 
fields for 77 = 1. Length scales of arrows and of x and y axes are arbitrary. More details can be found in |38 |. 


commutator terms which can be interpreted as effective 
chromo-electric and chromo-magnetic charges i?^], 

ig[A\,E\\ and ig[A\,B^, ig[A\,B\\, where El and B], 
are the transverse fields in nucleus k [5S] . The commuta¬ 
tors are non-zero when the gauge potential from nucleus 
1 can interact with the field of nucleus 2 and vice versa. 
Then for t > 0 the induced charges on opposite nuclei 
are indeed the negative of each other. Hence, we can 
also interpret the longitudinal fields as the abelian fields 
generated by additional color charges induced in the col¬ 
lision at t = 0. 

Going forward in time, we note that the first order 
in T brings no further contribution to the longitudinal 
fields, E^^ = 0 = but it is the leading order for the 
transverse fields 

Fw = -0 {[DloyFio)] ± ' (^ 6 ) 

Therefore, the transverse electric and magnetic fields 
grow linearly from their zero value at r = 0. We can 
express them in terms of the initial longitudinal fields as 


m 

M) =-I {smhg[D\Eo]+coshge^^[D^,Bo]) , 

(37) 

^( 1 ) = b (cosh?7 e*-’[H-’, ilo] - sinh77[I?*, Hq]) . 

(38) 


Recall that we have agreed to the notation H* = 9* — 
* 3 ^ 1 ( 0 )- [22] we have discussed extensively how these 

transverse fields can be understood from the QCD ana¬ 
logues of Faraday’s and Ampere’s Law for the rapidity 
even parts and from Gauss’ Law for the rapidity-odd 
components. In particular, it is very natural to expect 
rapidity-odd transverse fields even in a boost invariant 
situation. In Fig. we show a typical example for the 
rapidity-even and rapidity-odd initial transverse fields in 
an abelian example (covariant derivatives are replaced 
by ordinary derivatives in ( [3^ and (38l). Notice how at 
mid-rapidity field lines are closing around existing lon¬ 
gitudinal flux tubes (dark or light colored regions) due 
to Ampere’s and Faraday’s Law, while away from mid¬ 
rapidity Gauss’ Law allows for transverse flux between 



















longitudinal flux tubes. 

The first correction to the initial value of the longitu¬ 
dinal fields appears at order and in our short notation 
is 


Ff2) = ^[F^[F^Fo]], 

(39) 

Bf,~^ = ^[D\[D\Bo]]. 

(40) 


There is no correction to the transverse fields at this or¬ 
der, = 0. 

From order on the results become somewhat un¬ 
wieldy. For this reason we present the expressions for 
orders and in Appendix However, there is no 
particular reason why one could not in principle go to 
higher orders in powers of r. Generally, the longitudi¬ 
nal fields have only contributions for even powers of r 
and the transverse fields have contributions only for odd 
powers of r. 

To summarize this section, we have provided explicit 
formulas for the initial gluon field to an accuracy 

=FL„c + 0{t^) , (41) 

F* + 0{t^) , (42) 

for the electric field, and similarly for the magnetic field. 

III. THE ENERGY-MOMENTUM TENSOR OE 
THE FIELD 

From the field strength tensor we can easily calculate 
the energy-momentum tensor of the field 


The other diagonal elements of the energy-momentum 
tensor are 

T^(o) = 7^(0) = ^0 = -TfS) ■ (46) 

Hence the structure of the energy-momentum tensor for 
T —>■ 0 is the same as that for a longitudinal field in 
classical electrodynamics. There is a maximum pressure 
anisotropy between the transverse and longitudinal di¬ 
rections. Despite being far from equilibrium we take 
the liberty to use the notations of longitudinal pressure 
Pl = and transverse pressures iPx,Py) = 

We will denote the average transverse pressure as px = 

iPx+Py)/2. 

The initial transverse pressure pr = eo is large com¬ 
pared to an equilibrated system. A free, relativistic gas 
with the same energy density would have a transverse 
pressure £o/3. We expect a comparably large flow of en¬ 
ergy due to gradients in the transverse pressure. The lon¬ 
gitudinal pressure pl = —PT = —So is equally large and 
negative. The negative sign is not surprising. Keeping 
in mind the abelian reinterpretation of the longitudinal 
field, we expect the opposite sign induced color charges 
on the nuclei to be attractive. Hence, the initial lon¬ 
gitudinal fields would like to decelerate the sources. In 
fact, this is the mechanism that removes kinetic energy 
from the nuclei and deposits it as field strength in the 
space-time region between them. Here we do not take 
into account this back reaction of the field on the sources 
since the nuclei, even at top RHIC energies, seem to stay 
ultra-relativistic all the time, as discussed before. 

The qualitative global behavior of the system then is 
seemingly easy to predict from the simple form of at 
T —0, 


yMi- = p , (43) 

For brevity we will often employ a notation where 
SU{Nc) indices are summed over implicitly unless said 
otherwise: AB = A-B- = 2 Tr(Ai?), a = 1,..., — 1. 

We will now provide the first few orders in t for all com¬ 
ponents of the energy-momentum tensor 




(44) 


as functions of the initial longitudinal fields Eq and Bq. 


A. Initial Energy Density and Pressure 


/eo 0 0 0 \ 

0 eo 0 0 

(0) 0 0 eo 0 

Vo 0 0 -eo/ 


(47) 


While the negative longitudinal pressure leads to the de¬ 
celeration of the colliding nuclei, the transverse pressure 
forces the system to expand in the transverse direction. 
This transverse expansion, driven by the pressure of the 
classical field, is expected to be larger than in an equi¬ 
librated relativistic gas HH [IS]. We will see that this 
intuitive picture, while mostly correct, has to have addi¬ 
tional features added to it since the energy-momentum 
tensor above does not have the full information about the 
classical fields which drive the dynamics. 


Only the diagonal elements of have finite values at 
T = 0. We define eo to be the initial value for the energy 
density 

^o=T(°°) = 

= (lA\,Ai][AlAi]) . (45) 


B. Onset of Transverse Flow 

At the next order, linear in r, the components and 
T^\ with * = 1,2, are the only ones to pick up contri¬ 
butions. They describe the flow of energy and longitudi¬ 
nal momentum into the transverse direction. Note that 
is the transverse component of the Poynting vector 
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FIG. 4. (color online) Energy flow in the x-y-p\&ne in an abelian analogue for the same random distribution of seed fields 
Al, A 2 chosen in Fig. The transverse Poynting vector T'’* (arrows) is drawn with the energy density £0 (shading) in the 
background. The panels show rapidity-even (left), rapidity-odd (center) and full flow at r; = 1 (right). Length scales of arrows 
and of X and y axes are arbitrary. 


S = E X B. Therefore, the transverse expansion expected 
from the qualitative arguments given above sets in lin¬ 
early in r. We have 


Tth=^'^{B,El^-E,BQ 

(48) 

= cosh 77 -b ^, 0 * sinh 77 , 


rpZi jA rpi ID Toi 

^(1) -- ^O^(i) - ^O^(i) 

(49) 


= -a* sinhry -|- -/3* cosh 77 . 


We note that we have two contributions to transverse 
flow. The first term is the flow driven by the gradient 
of the transverse pressure as we would expect from a 
hydrodynamic picture [31] 

a* = -VVo . (50) 

The second term involves the 2-vector 




(51) 


The derivation of these and the following expressions is 
made easier by using a set of SU{3) identities assem¬ 
bled in Appendix [^ These flow terms have first been 
discussed by some of us in Ref. [33] . 

The /3® defies the naive expectations from our earlier 
analysis of the initial diagonal energy-momentum tensor. 
It is profoundly related to the electric and magnetic fields 
underlying the energy-momentum tensor. More precisely, 
it emerges from the rapidity-odd transverse fields man¬ 
dated by Gauss’ Law. The /3® enhances flow from larger 
to smaller energy densities in some regions and quenches 
it in other regions. This can be seen in the example of 
random abelian fields in Fig. [^ This abelian analogue is 
particularly interesting here since the non-abelian terms 
in (48) vanish in the event-average as discussed in [35] . 


However, they will be important when the field is sam¬ 
pled event-by-event. 

The contribution of to the energy flow is odd in 
space-time rapidity 77 . We want to stress that its ex¬ 
istence does not violate boost-invariance. Obviously /?* 
will have a role to play when angular momentum and 
directed flow in the system are studied. 


C. Order r^: Corrections to Energy Density and 
Pressure; Longitudinal Flow 


At order the diagonal elements of receive their 
first corrections and all the previously vanishing compo¬ 
nents acquire their leading contributions. On the other 
hand, the transverse flow of energy and longitudinal mo¬ 
mentum are not affected. 


T“) = 0 = r(f). 


(52) 


The expressions for the energy density, the longitudinal 
flow of energy, and the flow of longitudinal momentum 
are 


^00 

J(2) 


o oil 

ZT' ZT’O I D IDO I J^l rpl I TDl TDl 

^ 0 ^( 2 ) + ^ 0 ^( 2 ) + 

- ^ (W* + < 5 ) - ^ V*/ 3 * sinh 277 

-b ^ i 5 cosh 277, ( 53 ) 

8 


=-cosh 2?7 H— (5 sinh 2r ], 

8 8 


(54) 
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rp33 

-'( 2 ) 


TP IT’S D dS I ^ ZTii Tpi 

-^ 0 ^( 2 ) “ ^ 0 ^( 2 ) + 

i(V*a* + (5)-^V*/3* sinh 2 r; 
4 8 

+ i<5 cosh 277 . 

8 


1 

2 




(55) 


We have used Eqs. (A191 and (A20) to simplify these ex¬ 
pressions. Besides the divergence of the transverse fields, 
q;® and /?*, we find a new field that appears in the expres¬ 
sions above, namely 


At order we have 


rpOO 

J(4) 


_ TP Tp3 I TD 11)3 I Tpi TPi I TDl TDl 

—^0^(4) + -t>0^(4) + ^(1)^(3) + -^(l)-^(3) 

_i_ ^ IT’S jp3 I ^ d3 d3 
+ 2 (2)^(2) + 2^(2)^(2) 

=P+ cosh 277 -h ^crsinh 277. 


(63) 


<5=[Zl\Ao][l?*,Ao] + [Z?*,Bo][Z?\i3o]. (56) 

The divergence of the transverse flow is the expected re¬ 
action of the energy density to the initial flow, leading to 
depletion at the source and accumulation at the sink of 
the flow field. 

The remaining new contributions to this order give cor¬ 
rections to the transverse pressures 

rg,=i4^(E,\,E;., + B,\,B,\, (57) 

p2 tp2 d2 d2 a _i_ zt’ Z7’3 i td t:>3 
-il/(l)i^(l) - j + ^0^(2) + ^0^(2) 

= -i(-A,„ + S+(-l)‘u:), 

= (58) 

Here A is the 2-dimensional Laplace operator. There is 
no implicit summation over the double index i = 1 , 2 in 
the first equation. The new quantities are 

u;=^{[D\Eo]^-[D\Eo]^ (59) 

+ [D\Bo]^-[D^,Bo]^) , 

7 =[D\Eo][D^,Eo] + [D\Bo][D^,Bo ]. (60) 

The uj describes the anisotropy of the pressure in the x- 
and ^-directions and is therefore responsible for a phe¬ 
nomenon akin to elliptic flow in the transverse plane. 


7^03 _ ij 

J-ii) -e 


(■^a)^(3) + -^(3) 


WI) 


1 

32 


a cosh 277 -I- 


1 

32 


K sinh 277 , 


(64) 


7-133 _ 

^(4) - 


T7' 77’3 o o3 i i 07 

- -C/o£/( 4 ) - BqB^^^ + -I- 

1 


^ 77’3 77’3 ^ d3 o3 

- 2W2)W2) “ 2 (2)W2) 

— p -I- — K cosh 277 -I- — tj sinh 2rj , 


(65) 


T" ={-iy(El^El^)+Bl^Bl^^ 

tp2 jp2 d2 d2 a _i_ ZTi I7'3 i d td3 

“^(1)^(3) “ ^(1)^(3)J + ^0-^(4) + ^0^(4) 

_i_ ^ jp3 jp3 I ^ d3 d3 

+ 2 ( 2 )^( 2 ) + 2 (2) 

=P+{-iyX, (66) 


7-112 _ r^l 77’2 d1 d2 t^I t ^2 d1 t > 

-^[4) — ~ — n(^\r>, 


dl)W3) 


(1)^(3) 


d3)-^(l) 


d3)^(l) 
(67) 


D. Higher Orders 

At order the only contributions are the first correc¬ 
tions to the transverse flow and They are 


where the new coefficients p, k, a, A and v are explic¬ 
itly given in appendix]^ The expressions for the energy- 
momentum tensor discussed here are accurate up to cor¬ 
rections of order r® for the and components, and 
up to order r® for all other components. 


r“ =e*^- [BoEI^^ + Bf.^El^ - EoBf^^ - Ef.^B^ 
= ^ cosh 77 -I-C* sinh 77 ) , (61) 

T 3i TP jpi Tp 3 TPi TD TDl d3 JDl 

(3) - “ -^0^(3) - ^(2)^(1) “ ^0^(3) - ^(2)^(1) 

= —(^* sinh 77 -I-4 cosh 77 ) . (62) 
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We give the explicit expressions for the flow vectors 
and 4 iJ3 Appendix 


E. Checking Energy and Momentum Conservation 

The solutions of the Yang-Mills equations auto¬ 
matically satisfy energy and momentum conservation 
= 0. This can be checked explicitly order by or¬ 
der. The and receive contributions only 

for odd powers of r, whereas consists only of even 





powers. At order t we find, for ly = 0, 


IV. AVERAGING OVER COLOR SOURCES 
WITH TRANSVERSE DYNAMICS 


II 




+ 

= 0 , 


= (cosh„|--isinh,|)r” 

sink I, A + icoshijA^ r”, + VT“, 

-(V*a* +15) cosh ?7 

- cosh T] [—V*/3* sinh 2r] + 5 cosh 277 ] 

- sinh rj [—V*/3* cosh 2r] + 5 sinh 277 ] 

-V*Q!* cosh 77 + sinh 77 

( 68 ) 


and similarly for v = 3. 

Transverse momentum conservation, z/ = 1,2, is obvi¬ 
ous at zeroth order in r. From the corresponding equa¬ 
tion 




tO 


^cosh 77 — sinh 77 ^^ (69) 

- (^sinh77-cosh77|^)Tf*)+V*7)?*) 


=a* -b Wo , 


all terms containing the anomalous flow /3* drop out and 
the remaining expression obviously vanishes using the 
known result for the hydrodynamic flow a®. Note that 
the index i is not summed in the term containing T*®. 

At order we have a very similar picture 




= 3 cosh 77 — sinh 77 — 

' or] 


— 3 sinh 77 — cosh 77 — 

' or] 


rriSi 

J(3) 


rpOi 

^(3) 


+v®T(®) + 

= (-^£0 + -5 + (-l)®cc) -f V^ 7 ] ,(70) 


with the third order flow contribution C* dropping out. 
Again, the index i = 1,2 is not summed upon multiple 
appearance and in addition we define j to be the trans¬ 
verse index with j ^ i. Momentum conservation holds if 
the equation 


e = V®(-A£o + W(-l)®o7)-VJ7, (71) 

is true. It is proven explicitly in Appendix Similarly, 
the momentum conservation equations at order are 

64p + 3 k + V®f = 0 , 

3cr -b V®C* = 0 . (72) 

We are now confident that we have the correct analytic 
expressions for the initial gluon field. 


So far we have held the charge distributions pk in the 
two nuclei fixed. We have expressed the gluon fields and 
energy-momentum tensor after the collision in terms of 
the initial longitudinal gluon fields Eq and Bq and the 
initial transverse gauge potential Those, in turn, 

are determined by the gauge fields A^[pi] and ^ 2 (^ 2 ] in 
the two nuclei before the collision. In a given nuclear 
collision the color charge densities pk are not known to 
us. But if we know the statistical distribution of the 
densities we could use the results of the last two sec¬ 
tions for an event-by-event analysis in which color charges 
Pk are statistically sampled according to their distribu¬ 
tions. Averages over event samples can then be com¬ 
pared to event averages of experimental data taken. A 
CGC event generator of this kind, albeit in 2-blD, has re¬ 
cently been presented in the IP-Glasma framework ng. 
In that work the time evolution of the gluon fields in 
the forward light cone was solved numerically. An event 
generator based on our results would be able to sample 
fields or the energy-momentum tensor at early times di¬ 
rectly without solving differential equations. However, 
in this work we will rather focus on obtaining analytic 
results for the event averaged energy-momentum tensor. 
We use the assumptions of the MV model which pos¬ 
tulates a simple Gaussian distribution of color charges 
mis]. We have to generalize the MV model by allowing 
slowly varying average charge densities in the transverse 
plane. This will allow us to treat transverse gradients in 
pressure and their consequences. 

We start by observing that the expectation value of the 
color charge of any nucleus at any given point has to van¬ 
ish, {p{xj_)) = 0. However, we expect local fluctuations 
to occur on typical non-perturbative time scales which 
are much larger than the nuclear collision time. Hence, 
the fluctuations are frozen at the moment of the collision. 
The size of the fluctuations are given by the expectation 
value ]i ^ {p'^{x±)) of the squared charge density. In 
the MV model it is assumed that fluctuations are Gaus¬ 
sian, uncorrelated in space, and isotropic in SU{3). We 
will see later that it is necessary to introduce a finite 
resolution in space to regularize the UV divergence that 
would emerge from an infinite spatial resolution. When¬ 
ever taking averages (...) we will thus keep in mind that 
they have to be taken at a finite resolution. For an ob¬ 
servable O measured after the collision of two nuclei the 
expectation value is given by 


{^)pi,P2 


d[pi]d[p2]0{pi, P 2 )w{pi)w{p 2 ) , (73) 


where the weight functions w are Gaussians with widths 
given by the average local charge densities squared, pi 
and p 2 - 
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A. The MV Model with Transverse Gradients 

We start with a brief review of the MV model. We 
implement the averaging over color sources in a given 
nucleus by fixing the expectation values 

{pa(.x^,XA_)pb{y^.yiy) = ^ 

X A(a:'F,f_L)^(a;'F - ^ 

as a precise definition of a (light cone) volume density 
A(a::^,af_L) of sources for a nucleus moving along the + 
or — light cone. In addition, expectation values of any 
odd number of p-fields in this nucleus vanish. We have 
dropped the index k labeling a particular nucleus here for 
ease of notation, and a, b are explicit SU{3) indices. We 
have also made explicit the coupling constant g that was 
contained in p as defined in Eqs. (§ and The A (and 
p) are then volume (and area) number densities of color 
charge, summed over color degrees of freedom. Note that 
the normalization of A and p differ by a factor — 1 from 
many other occurrences in the literature, such as |39j . 
We allow for a dependence of the expectation value A on 
both the longitudinal coordinate and the transverse 
coordinate x±. 

The longitudinal smearing in x^ is necessary to com¬ 
pute expectation values correctly, as first realized in (TB] . 
A nucleus must be given a small, but finite, thickness 
across the light cone which we will do by introducing 

X{x^, x±) = p{x±)h{x^). (75) 

Here h is a, non-negative function with finite width 
around = 0 and normalized such that 

dx^A(x^,xj_) = p{x±). (76) 

It is not necessary to specify the shape of h further. 

We have introduced the dependence of the charge den¬ 
sities A and p on x±^ as a generalization of the original 
MV model, where the nuclei are assumed to be infinitely 
large in the transverse direction and on average invariant 
under rotations and translations. Real nuclei break these 
symmetries; in order to generate a non-trivial transverse 
dynamics we need to investigate how the results in the 
MV model generalize when small deviations from these 
symmetries are allowed. Our guiding principle is that, 
on transverse length scales that are equal to or smaller 
than the scale of color glass, l/Qs, the gluon field is de¬ 
scribed by the well-defined color glass formalism. On 
larger length scales other dynamical effects, for exam¬ 
ple from the nucleonic structure of the nucleus, appear 
and can be parameterized by the dependence of ^ on x_l . 
Here we introduce an infrared length scale l/m. We must 
require that p varies by a negligible amount on length 
scales smaller than l/m. Explicitly we require that 

\p{x±)\ > mT^ |VV(^_l)| 

» m-2 |V*VV(x_l)| > ... . (77) 


Then m is an infrared energy scale which separates color 
glass physics from long wavelength dynamics. It is nec¬ 
essary to have the hierarchy 

1/Qs < 1 /to <C , (78) 

where Ra is the nuclear radius. 

We have two main goals in this extended MV model: 
(i) Observables must be well behaved under small devia¬ 
tions from translational and rotational invariance, other¬ 
wise the original MV model would not be infrared safe. 
In practice this means that observables should be only 
weakly dependent on the infrared scale. We will explic¬ 
itly check this condition below, (ii) The results will allow 
simple long-wavelength dynamics, expressed in an expan¬ 
sion in gradients of p, which is compatible with color glass 
physics at small distances. In practice this will allow us 
to safely apply the MV model locally to realistic nuclei 
as long as the location is sufficiently far away from the 
surface of the nucleus where the density p starts to fall 
off quickly. 


B. The Gluon Distribution 

The most important expectation value of fields in a 
single nucleus is the two-point function {A{x±_)A{y _\_)) 
which, in light cone gauge, is related to the gluon distri¬ 
bution. The Yang-Mills equations ([^ for a single nucleus 
on the -f light cone are most easily solved in a covariant 
gauge first where A^^^ = The equations reduce to 

Aq;(x“,x_l) =-pcov(a;“,:?_L), (79) 

where the Laplace operator A acts on the transverse co¬ 
ordinates. The explicit solution is 

a{x~,xi_)= J dz'iG{xj_ - zj_)pcov{x~,zj_) , (80) 

with a Green’s function G{x±) = — ln(x^/A^)/47r where 
A is an arbitrary length scale. However, we will be bet¬ 
ter served by introducing a physically motivated regular¬ 
ization through a gluon mass m which can be inserted 
into the Fourier transformation of the Green’s function 
G{k) = 1/fc^ —>• l/{k'^ + w?) [35]. This gluon mass could 
be an unrelated infrared scale, but for simplicity we will 
choose it to be the same as the IR cutoff in the gradi¬ 
ent expansion of p introduced in the previous subsection. 
Including the gluon mass leads to the Green’s function 

G{x±) = ■^Ko{mx±) , (81) 

ZTT 

where Kq is a modified Bessel functions. This Green’s 
function reproduces the previous expression in the limit 
m —>■ 0 with A = 2e“'’'®/TO, where je is Euler’s constant. 
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The two-gluon correlation function in covariant gauge 
can then be easily derived from (74) as 


{aa{x ,xj_)ab{y ,yj_)} = 

X 6{x~ - y~)"/{x~,x±,y±) . (82) 

Here we have introduced another Green’s function 


l{x ,x±,y±) 

= j (PzxG{xx-zi.)G{yx- za_)\{x~,ZA_) ■ 

(83) 

We will see that 7 depends strongly on the IR regular¬ 
ization scale m. In the limit r = \y — x\ —>■ 0 it diverges 
like 1 /; cancellation of this divergence for observables 
is a critical test of the theory. 

The gluon field in light cone gauge can be derived 
from the gluon field in covariant gauge with the help of 
the Wilson line 


U{x ,x±) =V exp 


r 

—ig / a{z ,x±)dz 

J —00 


■ (84) 


Here V denotes path ordering of the helds a from right to 
left. One can show that the correct gauge transformation 
to arrive at the light cone gauge potential is [T5] 

AUx~, x±) = -U{x~,x±)d^UUx~,x±) . (85) 

9 

We apply this gauge transformation to the field strength 
tensors in covariant gauge to obtain the corresponding 
tensors in light cone gauge, F = UFcovU^- Their corre¬ 
lation function is 


This expectation value is suppressed since 'j{z~,x±,y±) 
tends to diverge in the limit m —> 0. For a double line 
we have 

{Ug^{x~,XA_)Ucd{x~^yx)) = 5gA5i^d{x~,xx,yA_), (89) 
where 


d{x ,f_L,yj_) = exp 


■ 

2{Ni - 1) 


X / dz r( 2 ; ,xx,yx) 
J —00 


is the exponentiation of the integral of 


(90) 


r(2: ,5 '_l,27j_) = 27(2 ,f_L,yj.) 

- 7 (z",f_L,f_L) - 7 ( 2 :",J/j_,y_L) (91) 

along the light cone. This F is a subtracted version of 7 . 
In the original MV model the subtraction in F removes 
the 1 /m^ singularity in 7 for small m and renders the 
exponential d finite. In particular, r(x“,x_L,y_L) van¬ 
ishes in the ultraviolet limit yx —)■ x_\_- We will show 
below that this crucial cancellation is still valid for our 
generalization. Here we have dropped contributions from 
non-color singlet pairs as in m- 

Now we return to the discussion of the correlation 
function of fields. One can prove that the only possi¬ 
ble contraction of fields on the right hand side of Eq. 

P comes from the factorization of expectation values 
'd^){d^ad^a) The second factor can be deter¬ 

mined from Eq. (821 as 


{d’^aaix ,x±)d^ab{y ,y±)) = 


9 




x6{x -y )VJ,V^ 7 (x ,x_l,27j_)- (92) 


Together with Eq. (89) this leads to the result 


{F+\x-,x^)F^\y-,yj_)) 

= (a;",f±) iy~,yx)) ■ 

( 86 ) 

In the above expression we have expressed the Wilson 
lines U by their counterparts in the adjoint representa¬ 
tion, U, by virtue of the relation 


{F+\x ,x±)F^^{y ,y±)) = _ ^ 5gx5{x -y ) 

X [V!,V^7(x",x_L,yj.)] <i(a;",x_L,27j_), (93) 

for the expectation value of fields in light cone gauge. The 
correlation function of two gauge potentials in light cone 
gauge follows from an integration with retarded bound¬ 
ary conditions 


UtaU^ = Uabh ■ 


(87) 


Let us take a small detour to discuss expectation values 
of adjoint, parallel Wilson lines in the MV model p^ . 
A systematic study was carried out by Eukushima and 
Hidaka [ID] . For a single line we obtain 


{Ugx{x ,f_L)) = ^Q^exp 


9^N, 


2 (V 2 - 1 ) 


/ 7(0 ,xx,x±)dz 

J —00 


■ ( 88 ) 


A*(x ,x_l) = — / dz F'^'^{z ,x_l). (94) 

J —00 

One integral is easily evaluated to give 
(AJx ,xi_)Al{y ,y_L)) = r ^ - 


X 


dx 


d 


dx'~ 


exp 


g^NJ'[xx,yx) 

9^N, 


2 (V 2 - 1 ) 


xr(x_L,y_L) / dz h{z ) 

J —00 


(95) 
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Note that we have used Eq. (75) which allows us to factor 
h{x~) from and r(xj_,y_L). We have formally 

defined ^{xx,y±) as the integral of 7 (a:“,xj_,yj.) over 
x~ from —oo to +oo, and similarly for E. Here we have 
rewritten one factor of h{x~) as a derivative dldx'~ of 
the exponential. 


We can now evaluate the second integral. We will 
only be interested in min{a;“,y“} > 0. Upon taking 
the limit of vanishing width of h we find that the fields 
are independent of the coordinates x~ and y~ as long as 
min{a;“,?/“} > 0. We simply write 


{A’;^{x±)Al{y±)) = 2g‘^5gh 




2x 


g^NJ'ixx^yx) 


X exp 


2 (iV2 - 1) 


r(f_L,y±) 


- 1 


(96) 


This result holds for both the MV model m and our 
generalization of it. 

Before proceeding, let us write down the correlation 
function of two gluon fields when we formally take the 
ultraviolet limit y± —> x±. In that limit T —>■ 0, and we 
can expand the exponential function around 0 , using only 
the two leading terms, to arrive at the simpler expression 


{Al{xx)Al{xx)) 


= S, 


ab 


9 


m-1 




y±^x± 


■ (97) 


For further evaluation of the gluon distribution we have 
to understand the correlation functions 7 and T. 


C. Gluon Fields in the MV Model with Transverse 
Gradients 


in Eq. 


to be 


7 o(?') = 7 o(^_l,? 7 j_) 

= Mo y d'^z±G{x± - zj_)G{yj_ - zj_) 


= Mo 


d k± _ 

(27r)2 (fc2_+m2)2 


/ 

Mo-j- Ki{mr ), 

47rm 


(98) 


where m is the same gluon mass introduced as a IR reg¬ 
ulator before. The 70 only depends on the relative dis¬ 
tance r = |r| due to isotropy and translational invari¬ 
ance. As mentioned before, 70 exhibits a quadratic de¬ 
pendence on the infrared cutoff m for small r, specifically 
it is 7o(0) = yolATTiTL^. 

On the other hand, this singularity cancels in the sub¬ 
tracted 2-point function (91). In the UV limit r —)■ 0 the 
leading term is 

ro(r) = 270 (r) - 270 ( 0 ) 


-I- 2"fE - 1 + Oirn^r"^). (99) 




This is the equivalent of the result in m using a finite 
gluon mass regularization. The Tq only exhibits a weak 
logarithmic dependence on m for small r. 

Let us now check that the same cancellation takes place 
if A is weakly varying on length scales 1/m as permitted 
here. We are only interested in typical values of r = 
\y± — ai_L| < <C m~^ since we will later take the UV 


limit. We recall from Eq. (81) that the Green functions 
G(z_l) ~ Ko{mz±) fall off on a scale l/m^ r. With this 
clear separation of length scales we can restrict ourselves 
to the first few terms of a Taylor expansion of fi around 
R in the calculation of 7 


y{zx) = y{R) + (zl - R)*VV(i?) + • • • • (100) 


This leads to 


The cancellation of the singularity in 7 through the 
subtraction in Eq. (91) is a classic result of the original 


MV model for constant (in transverse coordinates) aver¬ 
age charge densities. We will now show that this result 
holds for the inhomogeneous charge densities A that we 
have permitted. More precisely, we will show how expec¬ 
tation values of fields, like the gluon distribution above, 
can be systematically expanded in gradients of fi. Let us 
introduce center and relative coordinates for two points 
x± and y± in the transverse plane via R = {x± + y±)/2 
and r = y± — x±. The discussion in this subsection will 
use the area charge density g, but all results apply in 
a straightforward way to the generalized density A and 
correlation functions not integrated over x~. 

In the original MV model with constant y{x±) = /ig, 
we can easily calculate the correlation function 7 defined 


7(R,f) = 7(x_L,yj_) 

= lo{R,r) + ^V*VV(i?)7'' (7) + ... . (101) 

Here we have 70 (i?,r-) = fi{R)rKi{Tnr)/4:TTm analogous 
to Eq. ( [ 9 ^ , representing the constant term. The linear 
term vanishes because 

J d^zxG{zx + r/2)G{zx - r/2)z\ = 0 . (102) 

The second order term is 


1 "^ = J d'^zj_G{zj_)G{r- zj_)z]_z]_ 

= - :^K2{mr)-\ - —— - 

247 rm^ ASirm 


Ki{mr). (103) 
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These correlations functions can be conveniently com¬ 
puted in Fourier space, similar to the technique in Eq. 

([M). 

The subtraction of 7 ( 0 ) removes the leading quadrat- 
ically divergent term in m as in the original MV model. 
We can expand 70 and 7 *-^ for small mr. For F this leads 
to 


r(.R, f) = fi(R)— (inm^r^ ~ 2) -|- Oi^nn^r^) 

Stt ' ' 




-5^^ -f 


dSTrrn^ 


^ )+0(VV), (104) 


where m = mexp(jE + l/2)/2 « 1.47m. Indeed, the 
dependence on the cutoff m is at most logarithmic for 
the small variations of /j, that are permitted. Even though 
we could take the expansion (1011 farther we will never 


keep gradients of /r larger than second order. Higher 
derivatives will be hard to control phenomenologically, 
and it is now obvious that condition guarantees that 
the derivative correction in our result for F is small. 

Besides the subtracted correlation function F we need 
the double derivative Vl.V^yj(x, y) for the gluon distri¬ 
bution (97). As discussed above, we neglect gradients of 


/i beyond second order. We have two mass scales in the 
problem, m Q, which could cancel the dimensions of 
energy”^ introduced by the gradient expansion. The UV 
cutoff Q was introduced earlier as the resolution scale in 
the transverse plane. We anticipate that in the next step 
of the calculation we take the limit r —>■ 0 , meaning that 
explicit factors of r will turn into powers oil/Q. We only 
keep terms like m“^V ^ ^ rV. In other words, 

we drop terms that are suppressed by additional powers 
of the large scale Q. Thus we arrive at 


function with two derivatives in the ultraviolet limit are 


= —iJ.{R)—5'‘^ In(m^r^) 
Stt 


1 


dSTrm^ - 


r—>-0 

2V*VV(i?) -h Ay(R)5^^ 


(106) 


Equations (1041 and (106 1 , together with Eq. (96) with¬ 
out the gradient corrections, reproduce the standard re¬ 
sult for the 2-point function in the MV model mm 


{KixMAVE)) 


- 1 ) 
g2N^r‘^ 


1 - (m^r^)i6"(«?-i) 



(107) 


Remember that our definition of y has an additional fac¬ 
tor — 1 compared to Refs. |IH1|39]. 

Here we are strictly interested in the UV limit r —)■ 0 
regul arized by a resolution length scale l/Q. Plugging 
(106 1 directly into (97) we obtain 


(A^(£_L)A|(f_L)) = SgM 

v'=vV(^) 

QmP-yiix) 


gV(^) 

8 ^(V 2 - 1 ) 


02 

5^1 In ^ 


(^5^1 ^ij 25 *'= ( 5 ^') 


(108) 


keeping all leading terms in powers of l/Q up to sec¬ 
ond order in gradients. We have made the replacement 
r —>■ 1/Q in the logarithm, which is equivalent to impos¬ 
ing Q as the momentum cutoff in a Eourier representa¬ 
tion. The typical transverse momentum of gluons in the 
nuclear wave function is given by the saturation scale 
Qs- Here we can take Ql ~ g‘^y./{N^ — 1) in accordance 
with [35] (accounting for the factor — 1 difference in 
the definition of y). Qs is the largest scale in the prob¬ 
lem and thus the ultraviolet scale Q for a single nucleus 
should be proportional to Qs with some numerical factor, 

Q = KQs. 




= S'^^KQ(mr) — 

-*-* • • TT? ^ 

2V*VV(R)-f A/i(R)5*^ - :,KRmr) 

J 487 rm^ 


where any gradients V* on the right hand side act only 
on fJ.{R). Note that terms with single derivatives are 
power suppressed. Now we take the formal limit r —>■ 0. 
No dependence on the direction of r should remain in this 
limit and we keep only terms isotropic in f by setting 
—>■ 5*'l/2, The leading terms of the correlation 


D. Higher Twist Gluon Correlation Fnnctions 

For the components of the energy-momentum tensor 
beyond the leading term in the r-expansion, we will need 
expectation values of gluon fields beyond the 2 -point 
function. We will compute those correlation functions 
in this subsection. With more fields or more derivatives 
these are akin to higher twist distributions of the gluon 
field. The power counting technique in 1/Q we intro¬ 
duced in the previous subsection will be useful for book 
keeping. 

One additional transverse covariant derivative in the 
2 -gluon correlation function can be computed as follows. 
First, we again express gauge potentials in terms of field 
strengths 
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{D^A\{x ,XA_)Al{y ,yj_)) = f dx' f dy' /(£)'= F+*)a(x' ,XA_)F^\y' ,yx)) ■ 

J —oo J — oo ' ' 


(109) 


I 

Using the same change to covariant gauge as in Sec. tion value on the right hand side can be transformed into 


IV B and recalling that the expecta- the expression 


,y±)d''d"aa>ix' ,x±)d^ab'{y' ,y±) 


~y' ^x±,y±)]dix' ,x±,y±), ( 110 ) 


r 


in analogy to Eq. (93). Note that correlators with three straightforward to calculate. We have 


gluon fields vanish since an even number of adjoint Wil¬ 
son lines and fields a have to be contracted with each 
other. Combinations (Ua) ~ 0 are suppressed [lU] . 

The two integrals over x'~ and y~ can be dealt with 
exactly as in the case of the simple 2-point function. The 
result for arbitrary longitudinal positions a;“ > 0 (after 
taking the thickness of light cone sources to zero) is 


{D'^A^,{xx)Ai{xx)) 


m - 1 - 


(5^VLV^V^,7(^u,yu) 




( 111 ) 


in the interesting UV limit r —>■ 0. The same expectation 
value with the covariant derivative acting on the second 
gauge field would result in the same expression with the 
obvious replacement —>■ V^. 

We apply the same basic strategy to calculate expres¬ 
sions with more derivatives. We obtain 


{D’^A:{xx)D^AI{xx)) 


- 1 ) 


(112) 


In the same spirit we have 

{D'^D^Al{xx)Ai{xx)) 
g'^5gb 


- V^V^VlV^,7(^u,ja)- (113) 


The higher derivatives of the correlation function 7 are 


vLv^v^,7(^^,yu) = 


KR) 

47r 


+5^^—+ 5^^- ) mN:i(mr) 




-m^rK2{rnr) 


Stt V ^ / 


VV(i?) 

Stt 


_ ^Ugjk _ gkigij'^ Ko{mr), 


(114) 


where we kept the two leading orders, l/r and m, in our 
power counting in mr. One can check that the contri¬ 
bution of the leading term to observables, such as /3®, 
vanishes due to the odd number of powers in r*. Hence 
the relevant term in the UV limit is 




y±^x^ 


IGtt 


X In ± -f . (115) 

The lower signs are valid if the derivative V^' acts 
on y± instead of x±. The lower signs in the previ¬ 
ous expression will be useful for the expectation value 
{A'‘^D^Al{x±)). As a consistency check, we note that 


Eq. (115) switches between upper and lower signs under 


the exchange {i,a} O {j, 5} as dictated by symmetry. 
As discussed above, we have dropped a term 0{g‘^fiQ) 
that does not contribute to observables. 

Caution is needed when calculating four deriva¬ 
tives acting on 7 . The leading behavior of 

is similar to Alnr which 
vanishes everywhere except for r —>■ 0. A proper integra¬ 
tion will give us the leading term (again regularizing l/r 
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by Q) as 

= . (116) 

In the UV limit the next to leading term in the transverse 
scale hierarchy is 


'Q^ 




X In 


^2n 

^gij^km^ln _ ^ik^jmgln ^ gilpm^kn 


_|_ ^jkgimgln gjlgimgkn ^ ^kl^imgjn^ 


V-VXi?) fQ^\ 

967r \m'^ ) 

X ((jb jfcmji" _j_ gik^jmgln _|_ gilgjm^kn 
_|_ gjkgimgln _|_ gjl^imgkn _|_ gkl^imgjn'^ 

T — + 5^'^5^^ + 5^^5^^). (117) 

967r 


The expectation values of gluon fields with up to two 
covariant derivatives will enable us to calculate the ex¬ 
pectation values of components of the energy-momentum 
tensor up to order in the next section, including the 
effects of transverse flow. In addition, we will calculate 
energy density and pressure up to order To that end 
we also compute the leading Q'^ terms of the fourth order 
coefficients. However, we will neglect all effects of trans¬ 
verse gradients at fourth order, which would lead to very 
lengthy expressions. 


V. THE ENERGY-MOMENTUM TENSOR OF 
COLLIDING NUCLEI 


Note that in this chapter we calculate only averages of 
components of the energy-momentum tensor and will 
henceforth suppress the symbol (...) in the notation for 
simplicity. Applying (1081 for each nucleus, the initial 
energy density is 


2TTN^al (Qi\^ 

eo(a;_L) = I ^ 1 In 


Ql 

(119) 


where p,i and /i 2 are the expectation values of the densi¬ 
ties of charges in nuclei 1 and 2, respectively, and Qi and 
Q 2 are UV scales chosen for the wave function of nucleus 
1 and 2, respectively. We have dropped terms propor¬ 
tional to VV^/m^ which are subleading for the energy 
density. 


Expression (119) is very interesting. The appearance of 


as to the power 3 can be understood in the following way. 
This classical calculation corresponds to the emission of 
a gluon from source 1, the emission of another gluon from 
source 2, followed by their fusion via a triple gluon vertex. 
This involves 3 powers of the coupling g in the amplitude, 
hence to a power of 3 in when the amplitude is squared 
to get the energy density. The initial energy density is 
very sensitive to the numerical value of as, since changing 
it by a factor of 2 results in a change in the initial energy 
density of a factor of 8. Quantum corrections to the 
classical CGC results are difhcult to compute [221113] but 
may change this sensitivity dramatically. For example, 
it is reasonable to expect that one coupling is evaluated 
at the scale Qi, the second coupling at the scale Q 2 , and 
the third at a common scale Q. Using the lowest order 
renormalization group result for the running coupling 


as(Af') = 


1 


/32ln(MVA|cD) 


( 120 ) 


with (32 = (llA^c — 2A^/)/127r we would get 


After the discussion of gluon correlation functions in 
single nuclei we now return to the case of two colliding 
nuclei. We will further break down the expressions for 
the components of the energy-momentum tensor in the 
small T expansion in terms of the fields A\ and A 2 in the 
individual nuclei. It is then straightforward to apply the 
results of the last section. 


A. Energy Density and Flow 

The expectation value of the initial energy density Eq 
from Eq. (45) can be written as [39] 


^0 — (^ 0 ) 2 fabefcde 

(118) 


27rAcas(Q^) ^ ^ 

eo{x±)^ tJ-i{x±)Mx±). ( 121 ) 

This triumvirate of running couplings is reminiscent of 
what happens when computing quantum corrections to 
the small-a: evolution of the gluon distribution [44]. It 
appears that scale dependences are weaker once quan¬ 
tum corrections are established. Of course, the functions 
gi{x±) also depend to some degree on the scales. 

For phenomenological purposes we will introduce a 
common UV scale Q and, in what follows, we will al¬ 
ways make the simplification log((5f/m^) —>■ log(Q^/TO^), 
i = 1,2. For example, we can choose to be sim¬ 
ply the arithmetic mean of the scales of both nuclei, 
= {Ql -I- Q\)I2. This is a very good approxima¬ 
tion in the traditional MV setup where nuclei are con¬ 
sidered homogeneous slabs of color charges. For most re¬ 
alistic applications this will still be a reasonable choice. 
Recall that Qi is proportional to the saturation scale, 
Qi = KQsi for a given nucleus i, with a numerical factor 
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K. In collisions of two nuclei the relevant scale for the 
energy density is typically the larger of the two satura¬ 
tion scales However, experimentally accessible 

saturation scales do not cover a large range. Even for 
the largest nuclei at LHC energies they are at most a few 
GeV, barely one order of magnitude larger than Aqcd- 
Hence, assuming one common scale from some averaging 
procedure between both nuclei seems sufficient for many 
purposes. Because of the limited range in Qg we will also 
neglect a dependence of Q on the transverse coordinate 
which is, in principle, present. Thus we will not evaluate 
any transverse derivatives acting on Q. 

The expectation value of the rapidity-even flow vector 
in the transverse direction at order r is simply given by 


result 

HA\jD-^AAi,A,AiD^A)iAA')p2) 

+g^f^fcdeff^fh^ (< 5 *^< 5 '=' + 

MauAJp^ (AiA2,,AA^ATj}p2) ■ ( 128 ) 


a* = —27ra 


N, . , 


— —So 


" — 1 m 

V (MiA^ 2 ) 

AtlM 2 


In" (/X 1 M 2 ) 


( 122 ) 

(123) 


Separation of contributions from both nuclei for the 
rapidity-even flow vector leads to 


/?* = g^fg^fcdee^^ (124) 

X '{{D^AAAAlAiAhAj (125) 

+ {iAT.a)At){iD^AlAA‘)] ■ (126) 


The expectation value then takes a form complementary 

to a® [38] 


= -eo 


M 2 V >1 - MiV >2 

gig2 


(127) 


Note that the expectation value of /3* disappears for ^1 = 
/i 2 - Thus it vanishes for collisions of identical nuclei with 
impact parameter 6 = 0. We have discussed in detail in 
[35] how /3* describes a rotation of the fireball for 6 A 
0 while still preserving boost-invariance. We will come 
back to this in the next section. 


B. Higher Orders in r 


Using the higher twist gluon correlation function we de¬ 
rived in Sec. lIVDI this evaluates to 


S = 4£oQ^ In ^ 
+£o 




m‘ 


QA 


A^i A/r 2 

Ml m 2 


-|-147ra 




m - 1 ) 

147ra^ A^c 
' 3m2 (iV 2 _ 1 ) 

MiA ^2 + M 2 AM 1 
M 1 M 2 


^oln ( ^ ) (mi + M 2 ) 




eo 


2 (Ami + A^ 2 ) 


(129) 


The other coefficients at order are 


w 


7 


^^[v1v1(miM2)-V2v2(^i^2) 

-|-2 (V^Mi^^M2 ~ V^mi^^M2 )] 

Nc 57rag£o 
~ Nl-\ 3m2 

X [ - V1v1(mi + M 2 ) + V2v2(^i + /i 2 )] , (130) 

[Viv 2 (;ri^ 2 ) + (VV1VV2 + VV1VV2)] 

2 miM 2 


-k 


iV, 


107ras£o 

3to^ 


V^V^(Mi + M 2 ) ■ 


(131) 


It is interesting to note the hierarchy for terms at order 
of T^. Terms with derivatives are subleading to terms 
without, for example {tQA ^ (''■Q)^(VVm)/(to^m) ^ 
r^(VVM)/Mi while true non-abelian terms of order 
{rasAfi could be large as well. 

The energy flow A nt order can be expressed with 
the help of Eq. ( |7l] ) as derivatives of second order quan¬ 
tities. The leading correction to rapidity-odd flow at 
order is 


The expectation values of terms at order can be 
calculated in a straight forward but increasingly lengthy 
manner. Eor the coefficient S we have the intermediate 


c* 


^£ 0 ( 5 ^ In ^ 



MiV>2 - M2V>i 
M 1 M 2 

(132) 
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up to second order in transverse gradients. 

At fourth order in r we focus on the leading con¬ 
tributions for simplicity. For the relevant coefficients we 
obtain the expectation values 


3 ^4, -1 (Q^\ ^ 3 4 _2 /Q2 

P = In ^ In ^ 

32 J 8 \m‘‘ 

64 

« = -yp. 


(133) 

(134) 


VI. PHENOMENOLOGY OE CLASSICAL 
FIELDS IN HEAVY ION COLLISIONS 


With the results from the last section we are now ready 
to discuss the early time evolution of key quantities in 
high energy nuclear collisions analytically. We can com¬ 
pare to some numerical results available in the literature. 


A. Time Evolution of Energy Density and Pressure 


Let us begin by first considering the very simple case 
of homogeneous, equally thick nuclei, in other words, the 
case of colliding slabs with /j,i = /i 2 being constants. In 
that case any dynamics comes solely from the longitudi¬ 
nal expansion of the system. Because of its simplicity, 
this is an approximation often employed in the literature 
to study the general behavior of color glass systems. 

Neglecting transverse gradients, and keeping only the 
leading {tQ)^ terms, the results from the last section 
imply 


e = T°°(t,?7) = So 

, 3(Qr)4 
32a2 

+ 0{t^) 


1 - 


{Qry 


^1 — ^ cosh 277 ^ 


(a -I- 4) ~ ^ cosh 2 ? 7 ^ 


(135) 


PT = T^\t,T]) = So 
3(Qr)4 


1 - 


iQ'r) 

a 

^ 6 ' 


32a2 

PL = T^^{T,-q) = -So 
3(Qr)4 


+ :: (a + 4)+0(r^) 

{Qr? 


1 - 


(136) 

^1 -I- ^ cosh2?7^ 


32a2 


(a -I- 4) ^ 


(137) 


where we have defined a = h\{Q^/rh?) for brevity. We 
have neglected terms of order {asa)‘^Q^^ for two reasons. 
First, we have not computed the corresponding terms for 
the 4th order in time so we cannot evaluate these terms 
consistently. The calculation is somewhat tedious and 


reserved for a future publication. Secondly, the exact re¬ 
lation between and Q is not fixed from first principles. 
We can estimate that with a reasonable value K = 2 
the corrections to pr/s and Pl!^ are small up to about 
Qt « 0.8, which leads us to believe that the following 
analysis is valid. 

We can write down very simple but powerful pocket 
formulas for the time evolution of key quantities. For 
example the transverse and longitudinal pressure relative 
to the energy density at midrapidity behave as 




(138) 


Suppose we drop the order terms in the numerator 
and denominator. Then at a time given by (Qt)^ = |a 
we have pt = Pl = This corresponds to the equation 
of state of a massless gas of quarks and gluo ns. 

We can compare the results in Eq. (138) with those 
of Gelis and Epelbaum j^D]. They performed a real¬ 
time lattice simulation for colliding slabs using the gauge 
group SU{2). In Fig. [^we show results for the transverse 
and longitudinal pressures over the energy density, pr/s 
and Ph!^-, from our analytic approach up to fourth or¬ 
der in r and the numerical results from [20] (labeled LO 
in their work). Here we have chosen the values of a = 
0.8, 0.9, and 1.0, all of which give very good matching for 
small time, and are not unreasonable for small saturation 
scales. Note that this is a very schematic comparison for 
several reasons. A more quantitative statement would 
require a careful analysis of the IR and UV scales in the 
numerical calculation and their relation to Qs, a further 
investigation of p?" and terms in the analytic result. 



Qt 


FIG. 5. (color online) Evolution of the ratios of the transverse 
(upper curves) and longitudinal (lower curves) pressure over 
energy density for the classical gluon field at fourth order 
accuracy in time, compared to the numerical result from |20| 
at leading order for g — 0.5 (solid lines). Values of a = 0.8, 
0.9, and 1.0 (dash-dotted, dotted and dashed lines resp.) are 
indicated. 
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FIG. 6. (color online) Different ffow fields (black arrows) and initial energy density eo (shading) for Pb+Pb collisions at impact 
parameter 6 = 6 fm in the x — y-plane. The nucleus centered at a: = 3 fm travels in the positive ? 7 -direction. Left panel: a*. 
Center Panel: /?*. Right Panel: Full transverse Poynting vector at ?7 = 1. Note that a* is proportional to at = 0. 


and the use of SU{2) instead of SU{3) in our calcula¬ 
tions. However, it is interesting to note that the results 
agree quite well up to r ^ 1/Q. 

The comparison with numerical work is important in 
two ways. First, the study in [20] indicates that clas¬ 
sical field dynamics is sufficient for times smaller than 
To ~ 1/Qs, at least at small to moderate values of the 
strong coupling g. After that time quantum corrections 
and instabilities start to dominate. The successful com¬ 
parison also validates our previous argument about the 
convergence radius of the small-time expansion which we 
expected to be given parametrically by l/Qg. Indeed we 
can reproduce the results for transverse and longitudinal 
pressure very well up to that time. If we would want to 
relax the conditions and allow transverse gradients, we 
would also introduce dimensionless terms rV* which are 
smaller than tQ in the region of applicability. 

Serendipitously our near-held expansion works rather 
well up to the same time scale to which the classical held 
approach is valid. Thus we are led to believe that our 
analytic results are a rather simple and almost complete 
account of the collision dynamics up to tq. The asymp¬ 
totic values for pt/s and pl/s reached in the classical 
theory after tq are ~ 1/2 and ^ 0, respectively. Quan¬ 
tum corrections and instabilities will, however, lead to 
further isotropization soon after tq [20 ] . 


B. Global Flow of Glasma 

Two of us have discussed the effect of the two hrst- 
order how terms a* and /I* in detail in [38] . The 
hydrodynamic-like how term a* obviously leads to both 
radial and elliptic how; see left panel of Fig. 6. Note that 
this is how of energy of the classical gluon held at this 
point. However, due to energy and momentum conserva¬ 
tion, this how will translate into a how of huid cells after 
thermalization. We will discuss this in a future publica¬ 


tion. 

The rapidity-odd how term /3* potentially has many 
interesting implications; see center panel of Fig. 6. Its 
event average vanishes for central collisions (impact pa¬ 
rameter 6 = 0) for collisions of identical nuclei. However 
for hnite impact parameters it carries the angular mo¬ 
mentum of the gluon held that is transferred from the 
non-vanishing angular momentum of the two colliding 
nuclei. The how held exhibits a characteristic rotation 
pattern around the impact vector; see the right panel 
of Fig. 6 and top panel of Fig. [^ This would lead 
to directed how vi of particles which has been observed 
in experiments. The angular momentum would be trans¬ 
ferred to the quark-gluon huid at a later stage with poten¬ 
tial interesting consequences [IZIIIHI. We again refer the 
reader to [38] for more details. In collisions of two differ¬ 
ent species of nuclei, /3* leads to an increase of the radial 
how in the wake of the larger nucleus while suppressing 
how in the wake of the smaller nucleus; see bottom panel 
of Fig. [^ For asymmetric collisions at hnite impact pa¬ 
rameter the how held becomes more complicated. This 
could lead to interesting how patterns unique to classical 
gluon held dynamics [3H|. Those could be a novel signa¬ 
ture for the importance of color glass condensate in this 
regime. For the illustrations shown here, Woods-Saxon 
prohles have been used for the volume density of nucle¬ 
ons in the nuclei from which the transverse color charge 
densities /ri ^2 are computed. 

The second order in time also introduces a pressure 
anisotropy in the transverse plane for asymmetric colli¬ 
sion systems. The eccentricity of the transverse pressure 
Cp = (T^^ — T^^)/(r^^ -I- is often used to measure 
the buildup of elliptic how in the system. For the event 
average we read off from Eq. (57) that 


(^pix,y,v) = 


uj{x,y)T'^ 


(139) 
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up to second order in gradients and up to second order 
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FIG. 7. (color online) Transverse Poynting vector (black 
arrows) and initial energy density eo (shading) in the rj — x 
plane at t/ = 0. Top panel: Pb+Pb at fo = 6 fm. Angu¬ 
lar momentum is carried by the gluon field. Bottom panel: 
Pb-|-Ca at & = 0 fm, the Pb nucleus is moving to the right. 
The system expands more strongly in the wake of the larger 
nucleus. 


in T. This quantity is independent of i]. We see that the 
pressure anisotropy indeed starts to grow quadratically 
in time. We leave further numerical analysis to a future 
paper. 

Third order corrections typically slow the linear growth 
of the energy flow. For the rapidity-even part we again 
have a compact formula if we neglect terms with three or 
more derivatives. From Eq. and the expression for 
the expectation values of a* and <5 we obtain 

TeLn = ^ ^ cosh T] . (140) 

Similarly, from the expectation values for 13'^ and C i 
have 

(l - (141) 



FIG. 8. (color online) Time evolution for \V\ (solid black line) 
and for the radial projection of T®* (arbitrary units, dashed 
red line) for central Pb-(-Pb collisions at midrapidity with 
approximations as described in the text. 
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FIG. 9. Dependence of \V\ on the radial coordinate r for 
central Pb-|-Pb collisions at midrapidity. The radius of a Pb 
nucleus is about 7 fm. 


when higher order gradients and terms or order /i? are 
neglected. Interestingly, when we look at E* = 
at midrapidity as a proxy for velocity, the leading cor¬ 
rections in the time evolution cancel in numerator and 
denominator. They are of order /2a for both the 

energy density and the even part of the energy flow. In 
other words, while the growth of slows and invariably 
peaks and diminishes due to the longitudinal expansion, 
the velocity E* continues to grow roughly linearly as 

+ 0{t^) , (142) 

2 £o 

at midrapidity when transverse gradients of third order 
and higher are neglected. Figureshows the time evolu¬ 
tion of the radial velocity \V\ up to corrections of order 
T® for a point r = 5 fm away from the center of a central 
Pb-l-Pb collision at midrapidity. We also computed the 
time evolution of the radial projection of including 
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the correction, to contrast its slowing down to the lin- of this calculation, 
ear growth of . For the calculation of we have 
chosen Q'^ = 2 GeV^ and a = 1. Figure displays the 

radial dependence of |y| for the same central Pb+Pb col- C. Towards Quark Gluon Plasma 

lisions at T = 0.2 fm. We see that the surface velocity 

peaks around 0.2. However, one has to be cautioned that Let us summarize our knowledge of nuclear collisions 
typically the first fermi of the boundary (beyond r = 6 at a typical time tq = l/Qs- The energy-momentum 
fm for a Pb nucleus) is usually outside of the applicability tensor can be written, up to third order in r as 


T 


mn 


/£o-^(-2Aeo + 5) + 

f -f £o - X (-Aeo + 5 - w) 

+ 7 


Here we have used the r, cc, y, y coordinate system for the 
tensor. This gets rid of unwieldy coshy and sinhy terms 
from boosts. Note that there is no explicit dependence 
on y in this coordinate system due to boost invariance. 
This tensor exhibits the standard features expected of a 
fireball: radial and elliptic flow, and a decrease of en¬ 
ergy density and pressure with time, mostly due to the 
longitudinal expansion. In addition, we find angular mo¬ 
mentum and directed flow for finite impact parameter 
collisions, and a complicated flow pattern for asymmet¬ 
ric collision systems. These features can be predicted 
more or less accurately and in analytic form averaged 
over events. 

The reader should keep in mind that the phenomeno¬ 
logical analyses in the present section are rather crude 
and could be refined in many ways, as pointed out nu¬ 
merous times. However, they result in compact pocket 
formulas which could be useful for quick estimates in 
many situations. A more careful analysis can be done 
starting with the full expressions from Sec. [V] 

After a proper time Tq, instabilities growing from small 
fluctuations take over, leading to turbulent behavior of 
the fields. Further isotropization and equilibration is 
then expected to lead to quark-gluon plasma near ki¬ 
netic equilibrium. From a phenomenological perspective, 
one could simply translate the energy-momentum tensor 
of the classical field around the time tq directly into hy¬ 
drodynamic fields, as was done in m for ideal hydrody¬ 
namics and in El for viscous hydrodynamics. However, 
this obviously leads to large shear stress corrections, as 
can be seen from the large difference between transverse 
and longitudinal pressure around tq as presented previ¬ 
ously. It would be very interesting to see how key features 
of the transverse flow field translate into hydrodynam¬ 
ics and how they fare during subsequent hydrodynamic 
evolution. This would enable us to connect features of 
classical gluon fields in the initial state to observables. 

It would be relatively straight forward to build a semi- 
analytic event generator from our results. For example, 
one could follow reference m which used a model for 






So ~ ^i~^so -|- (5 -l- w) 


5/3" 

Ipy 

|(-2A£o + 3<5)/ 


(143) 


charge configurations of nuclei in collisions. In our ap¬ 
proach, their numerical solution to the Yang-Mills equa¬ 
tions would be replaced by our analytic time evolution 
using the near-held approximation. Then, from the sam¬ 
pled charge distributions, one has to calculate the co¬ 
efficients Co, a®, /3®, 5, etc. to obtain an event-by-event 
energy-momentum tensor. 


VII. CONCLUSION 

In this paper we worked out analytic solutions of the 
Yang-Mills equations for two nuclei with random color 
charges colliding on the light cone. Using a recursive so¬ 
lution we computed the early time gluon held and energy- 
momentum tensor in a near-held approximation. We hnd 
that this approximation gives acceptable results roughly 
up to a time Tq given by the inverse of the saturation scale 
Qs ■ This coincides with the time at which the entire clas¬ 
sical held approximation starts to breaks down anyway. 
Explicit expressions for the helds and energy-momentum 
tensor up to order have been provided. 

We have also calculated expectation values for the 
energy-momentum tensor when many events are av¬ 
eraged. Our calculation generalizes the McLerran- 
Venugopalan model to allow small but non-vanishing gra¬ 
dients in the average color charge in the transverse plane. 
This permitted us to discuss how phenomena in aver¬ 
aged events. We provide a comprehensive set of expecta¬ 
tion values of coefficients of the energy-momentum tensor 
which allow predictions for event-averaged for times 
around tq. We give compact and analytic formulas for 
key quantities like the time evolution of energy density, 
transverse and longitudinal pressure, the time evolution 
of transverse how of energy, and the time evolution of the 
transverse pressure asymmetry. 

We hnd that the transverse how of energy grows lin¬ 
early with time and that it can reach sizeable values at 
the surface of the hreball at tq. We have also discov¬ 
ered that the asymmetry between transverse pressures 
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starts to grow quadratically in time. The time evolu¬ 
tion of transverse and longitudinal pressure matches well 
with numerical results available in the literature up to 
tq. Besides the usual radial and elliptic flow a rapidity- 
odd flow emerges. We suggest that this energy flow of 
the glasma could be the origin of directed flow. It car¬ 
ries angular momentum which rotates the fireball. More 
complex flow patterns appear for collisions of asymmet¬ 
ric nuclei. The characteristic glasma flow pattern could 
potentially lead to another signature for color glass dy¬ 
namics in high energy collisions. 

At To ~ 1/Qs our calculation becomes unreliable. 
However, it could be attempted to match our results to 
a (3-|-l)-D viscous hydrodynamic code. We will discuss 
this in a forthcoming publication. We have also discussed 
the possibility to construct an event generator based on 
the results of this paper. 


It is useful to express Cartesian and light cone derivatives 
via hyperbolic ones by 


± d 


T 


1 d 


T dr 2x^ drj 


(A4) 


and 


d , d I . , (9 

— =cosh?7--smhry— , 

ot OT T orj 

d , , 9 1 , 9 

— = — smh?7-— I — coshry— . 

oz OT T orj 


(A5) 

(A6) 


Our conventions for covariant derivatives and field 
strength tensors are 
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Df^ = df^- igA ^^, (A7) 

- d''A>^-ig[Af^, A’']. 

^ (A8) 

Here A^, and are SU{Nc) valued functions that 
can be expressed as linear combinations of the SU{Nc) 
generators t-, a = I,..., — I. The generators are 

defined through [t-, t-] = if-t- and normalized by 

Tr = 0 , ^• (A9) 

This immediately implies that 


Appendix A: General Definitions 

Some conventions and useful formulae are gathered in 
this appendix. 3-vectors are denoted by bold symbols, 
vector arrows denote 2-vectors in the transverse plane. 
As an example, = (t,x) = (t,x±,z). Light cone coor¬ 
dinates are defined by 

^ {x° ± x^) , (AI) 

with d^x = dx'^dx~(Px± and = x^y~ + x~y^ — 

x\y\. Note that = {d/dt,—V) and = djdx^. 
Unless indicated otherwise, small Latin indices f,j, A: in¬ 
dicate transverse components of a vector, Greek indices 
label 4-vectors in (t, x, y, z) coordinates, and Latin in¬ 
dices m,n label 4-vectors in {T,x,y,7j) coordinates. Un¬ 
derlined Latin indices refer to the SU{3) algebra. 

Proper time t and space-time rapidity rj for a space- 
time point x'^ are defined as 


T = \/F — z^ = \/2x'^x~ , 



Tr(Y)=0, Tr([Y,Y]) = 0, (AIO) 

for any Y, Y in the SU{Nc) algebra since [Y, Y] G 
SU{N,). 

Using the ordinary product rule and the Jacobi iden¬ 
tity, one can show that covariant derivatives obey the 
generalized product rule 

[D\XY] = [D^X]Y + X[D\Y], (All) 

for any Y, Y in the SU{Nc) algebra; in particular 

[D\ [Y, Y]] = [[D^X],Y] + [Y, [D\ Y]]. (A12) 

It is sometimes helpful to interpret Ti{AB) as a bi¬ 
linear scalar product on SU{3) and [A,B] as a skew- 
symmetric product whose result is orthogonal to both Y 
and Y such that 

Tr(Y,[Y,Y]) = 0. (A13) 


(A2) 

(A3) 


This leads to some important ways to simplify expres- 
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sions with a trace involved. They include: 


Tt:{X[D\X]) = 
e^j[D\[Di^X]] = 
Tt:{[D^X][D^X]) = 

Tr{^i[D\X][D\Y]) = 

X 

+ 

Tt{X[D\[D\X]]) = 


-iv*Tr(X2), (A14) 

ig[B\X], (A15) 

^ATr(A 2 ) 

1:t{X[D\[D\X]]) ,(A16) 
2 

Tr([i4^A]y- [£>^y]A) 
]^y^t:v{[d\[d\x]y 
[D\[D\Y]X), (A17) 

^ATr(A 2 ) 

Tt{[D\X][D\X]), (A18) 


where X and Y are any SU{Nc) fields. As an example, 
the second equation implies that 


e*^■Tr(^;o[44^[i^^i?o]]=0, (A19) 

e^^TT{Bo[D\[D^,Eo]]=0. (A20) 


Appendix B: Expressions at Order r® and 

At order the transverse fields are 

+ |(e*^[Bo,T^5±[Ao,F(\±]) 

T |0e*^■[74^[^;o,Bo]], (Bl) 


whereas = 0 = B^^y In terms of the initial fields 
the third order fields are 


Els) = - ^ (coshr,e*^■[Z4^ [Z4^ [Z4^i?o]]] 

+ sinh77[Z4^[I4^[i4^^;o]]]) 

- sinhr? {[Bo, [D^Eo]] + [Eo, [D^,Bo]]) 

- ^ coshr; {[Eo, D\ Eo] - [Bo, [D\ Bo]]) 

-^^smhgY^[D^,[Eo,Bo]], (B2) 


Blo^=-^{smhr,[D\[D\[D\Bo]]] 

-coshrjY^[D^,[D^[D'‘,Eo]]]) 

- ^ cosh 77 {[Bo, [D\ Ao]] + [Eo, [D\ Bq]]) 

- ^ sinhr 7 e*^- {[Bo, [D\Bo]] - [Eo, [D^,Eo]]) 

-'^coshr,[D\[Eo,Bo]]. (B3) 

The longitudinal field at order r'* is 
Ef,^=^^[D\[D\[D^,[D^,Eo]]]] 

+ ^ige^^[[D\Eo],[D^,Bo]], (B4) 


Bf,)=^[D\[D\[D^,[D^,Bo]]]] 

-^ige^n[D\Eo],[D^,Eo]] 

+ ^^gY^[[D\Bo],[D^,Bo]] 

+ ^[Eo,[Bo,Eo]]. (B5) 


For the energy-momentum tensor t he t ransverse flow 
vectors and C*, as defined in Eq. (621, are given in 
terms of Eo and Bo by 


e = [D\Eo[D^,[D\Eo]]+Bo[D\[D\Bo]]] 

+ [D\Eo][D‘', [D‘,Eo]] + [D\Bo][D\[D\Bo]] 

- igY^Bo[Eo,[D^,Eo]], (B6) 

C = {[D\Eo[D\[D^,Bo]]-Bo[D^,[D\Eo]]] 

- 3[D^, Eo] [D^ ,[D\Bo]]+ 3[D^ , Bo] [D‘,[D^,Eo]]) 

- 3tgEo[Bo,[D\Bo]]. (B7) 

The components which we defined at order are 


P —iIoB(4) + Aoi?(4) + -(B(2)B(2) -I- A(2)A(2)) , 


(B8) 


n=[D\Bo][D\[D\[D^Bo]]] (B9) 

+ [D\Eo][D\[D'^,[D^Eo]] 

+ igY^[D\ Bo]{[Bo, [D\Bo]] - [Eo, [D^,Eo]]) 
+ igY^[D\ Eo]{[Bo, [D\Eo]] + [Eo, [D\Bo]] 

+ [D\[Eo,Bo]]]), 


a =Yi[D\Eo][D\ [D\ [D\ Bo]]] (BIO) 

-Y^[D\BoW,[D'^,[D>^,Eo]]] 

+ ig[D\Bo]{[Bo,[D\Eo]] 

+ [Bq, [D’', Bq]] -|- [D^, [Eo, Bo]]) 

+ ig[D\ Eo]{[Eo, [D\ Eo]] - [Bo, [D\ Bo]]), 
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\ TT'l j^l I dI d1 Z7'2 77'2 D2 d2 

A + ^(1)'^(3) ~ '^(1)^(3) ~ ^(1)'^(3) ’ 


(Bll) 


,, _ pi p2 pi r 2 pi p2 pi p2 

— “^( 1 )^( 3 ) “ -°( 1 )'°( 3 ) “ '^( 3 )^( 1 ) “ -°( 3 )'°( 1 ) 


(B12) 


We omit the lengthy expression for A and v in terms of 
Eq and Bq. 


Appendix C: Energy-Momentum Conservation at 
Order r® 


the proof for would be similar. 

Vi(-Aeo + i5-w)- 

= - V\{{Eo[D\ [D\Eo]] + Bo[D^, [D\Bo]]) 

+ i[D\ Eo] [D ^, Eo] + [D ‘, Bo] [D ‘, Bq]) 

- ([B',Bo][B',Bo] + [B',Bo][B',Bo])) 

+ ^{[D\[D\Eo]^] + [D\[D\Bo]^] 
-[DMD^,Eof]-[D\[D^,Bo]^]) 

+ [D^, [D\Eo][D^,Eo] + [B\Bo][B2,Bo]] 

= [B\ Bo [B', [D ‘, Boj] + Bo [D ^, [B' , Bq]]] 

+ [B\ Bo] [B', [B', Bo]] + [Bi, Bo] [B*, [B', Bo]] 

• (Cl) 


We prove explicitly the conservation of tr ansverse mo¬ 
mentum at order r^, i.e. the equation ( (tT] ). We do this 
for the first component (—Aeo + <5 — w) — V^y, 


Here we have used the product rule for covariant deriva¬ 
tives ext ensive ly and Eq. (A16) for the first equal sign 
and Eq. (A15) at the second equal sign. 
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